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Abstract. - Recently [EPL, 86, (2009) 17001], we had simulated the classical Langevin dynamics 
of a charged particle on the surface of a sphere in the presence of an externally applied magnetic 
field, and found a finite value for the orbital diamagnetic moment in the long-time limit. This 
result is surprising in that it seems to violate the classic Bohr-van Leeuwen Theorem on the 
absence of classical diamanetism. It was indeed questioned by some workers [EPL, 89, (2010) 
37001] who verified that the Fokker-Planck (FP) equation derived from our Langevin equation, 
was satisfied by the classical canonical density in the steady state, obtained by setting d/dt = 
in the FP equation. Inasmuch as the canonical density does not contain the magnetic field, they 
concluded that the diamagnetic moment must be zero. The purpose of this note is to show that 
this argument and the conclusion are invalid - instead of setting d/dt = 0, one must first obtain 
the fundamental time-dependent solution for the FP equation, and then calculate the expectation 
value of the diamagnetic moment, and finally consider its long-time limit (i.e., t — >■ oo). This 
would indeed correspond to our numerical simulation of the dynamics. That this is indeed so is 
shown by considering the simpler analytically solvable problem, namely that for an unbounded 
plane for which the above procedure can be carried out exactly. We then find that the limiting 
value for the expectation of the diamagnetic moment is indeed non-zero, and yet the steady-state 
FP equation obtained by setting d/dt = is satisfied by the canonical density. Admittedly, the 
exact analytical solution for the sphere is not available. But, the exact solution obtained for the 
case of the unbounded 2D-plane illustrates our point all right. We also present some further new 
results for other finite but unbounded surfaces such the ellipsoids of revolution. 



Absence of classical orbital diamagnetism - the Bohr- 
van Leeuwen (BvL) theorem has been a surprise of 
physics [Sjl^ . This is so because the amperean current as- 
sociated with the closed classical orbits of a charged parti- 
cle under the influence of the Lorentz force is expected to 
subtend a non-zero moment that should be diamagnetic 
- the Lenz's law. Bohr had given a physical argument 
for this null result in terms of the cancellation of the dia- 
magnetic moment of the completed orbits of the charged 
particle away from the boundary by the paramagnetic mo- 
ment due to the incomplete orbits skipping the boundary 
in the opposite sense. This surprising cancellation had mo- 
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tivated us in a previous work [5j to examine the case of a 
charged particle moving on a finite but unbounded closed 
surface, namely, the surface of a sphere, that naturally has 
no boundary. To our surprise, the long-time average of the 
orbital motion indeed gave a finite diamagnetic moment. 
The motion was treated through the stochastic Langevin 
dynamics that is expected to give equilibrium in the long- 
time limit - the Einsteinian approach to statistical me- 
chanics 1 6 . The unbounded nature of the sphere ensured 
avoided cancellation. The above result obtained by us has 
been recently questioned by some workers |7 who verified 
that the Fokker-Planck (FP) equation derived from our 
Langevin dynamics was satisfied by the classical canonical 
density in the steady state obtained by setting d/dt = in 
the FP equation. Inasmuch as the canonical density does 
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not contain the magnetic field, they concluded that the 
diamagnetic moment is zero. The purpose of this note is 
to show that this conclusion is invalid - instead of setting 
d/dt = 0, one must first obtain the fundamental time- 
dependent solution for the FP equation, then calculate 
the expectation value of the diamagnetic moment, and fi- 
nally consider its long-time limit (i.e., t — oo) - this would 
indeed correspond to our numerical simulation of the dy- 
namics. That this is indeed so is shown by considering the 
simpler but, analytically solvable problem, namely that for 
an unbounded plane for which the above procedure can be 
carried out exactly. We then find that the limiting value 
for the expectation of the diamagnetic moment is indeed 
non-zero, and yet the steady-state FP equation obtained 
by setting d/dt = is satisfied by the canonical density. 
Admittedly, the exact analytical solution for the sphere is 
not available. But, the exact solution for the case of the 
unbounded 2D-plane illustrates our point all right. We 
also present some further new results for other finite but 
unbounded surfaces such as the ellipsoids of revolution. 

Let us consider a charged particle (of charge — |e| and 
mass m) moving with a velocity u = {ur^^Uy) in a two- 
dimensional unbounded plane in the presence of a uniform 
magnetic field B perpendicular to the plane and under 
the influence of a stochastic force and the concomittant 
dissipation as in the classical Langevin equation 



u = — 7U u X B + /(t), 

mc 



(1) 



where 7 is the frictional relaxation rate, f{t) is a Gaussian 
white noise satisfying 



(2) 



The corresponding Fokker-Planck (FP) equation can now 
be written as [Sj 

aP _ _ dP_ dP d[{ju^+uj,Uy)P] 
dt dx ^ dy dux 

d[{-iUy - i^cUx)P] -fksTfd^P d^P 
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for the normalized transition probability density 
P(r, u, t|ro, Uo) that describes the stochastic evolu- 
tion of the system, where r = {x^y) is the position vector 
of the particle. Following the method due originally to 
Chandrasekhar |9j, the exact analytical solution of the 
above FP equation can be readily written down [sj as 

P = AAexp |-^^ (^^lup -2Hrn 

+2i7— (ux r)2 + CiG|r|' 
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Wexpj-^^^ [f{uI + ul) ~ 2H{xUx + yuy) 
+2H'^{yUx~xUy) + CiG{x^ +y^))Y (4) 



where without loss of generality we have set the initial 
conditions r,, = Uq = 0, i.e., a particle starting from rest 
at the origin. In the above equation Uc — \e\B/mc is the 
cyclotron frequency and 
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Note that equation Q is not invariant under field reversal 
B -B. 

Using the above solution for the FP equation, we can 
readily calculate the expectation of the orbital magnetic 
moment 
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Note that the expression for (m) is odd in B (or ujc)- In 
the long-time limit 



lim {m(t)) 



lei ksTLOc 



t-¥oo' ~ " mc (7^+0;^) 



(12) 



(This is precisely the expression given in equation (8) of 
[To]). It may also be noted from equation (111 that the 



time-scale involved in the growth of the magnetic moment 
from zero to full value is order 7^^ which is a micro- 
scopic quantity. However, the diamagnetic moment gets 
distributed over the whole space on a much longer time- 
scale. 

We now return to the main point of our contention. 
First, following [t], we verify that the classical (canonical) 
Boltzmann equilibrium distribution namely 
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2kBT 



{^1 + ul) 



(13) 



when substituted on the right hand side of equation ^ 
does give zero. This corresponds to a statement made in 



p-2 



Classical diamagnetism 



[7] that the classical canonical distribution satisfies their 
FP equation with d/dt = 0, i.e., steady state. The exact 
time-dependent solution, however, leads to a non-zero dia- 
magnetic moment in the limit t — >■ cx) as seen from equa- 
tions pTj ) and ( 12 ). This clearly shows that the arguement 
of [7] leading to the conclusion that classical diamagnetism 
is zero, is invalid. 

To reemphasize this point further in the context of the 
motion on a sphere, we would like to note that even the 
Liouville equation of |7| for a particle moving on the sur- 
face of a sphere has other solutions too and with non-zero 
currents. In fact, using Mathematica 
(notation as in [?]) 



11 , we found that 



Pst = sin 6* J"(ci,C2), 



where ci and C2 are given by 



ci = sind — -aujccos^ 



(14) 



(15) 



C2 



- [4(8«,2 -,- + aW^) + 4i4vl - aW^) cos 29 



-a Wc(32ti0 sin 9 cos^ 



a uJc cos 



(16) 



and T an arbitrary Junction of its two arguments ci and 
C2, satisfies equation (8) of [?]. And yet the solutions, 
equation (14 1, too have a finite magnetic moment inas- 



much as they are not invariant under the inversion of field 
B ^ -B. 

The cases of unbounded systems studied above clearly 
point to the role of the boundary in the context of the 
BvL theorem in that the absence of the boundary gives 
the avoided cancellation, and hence the non-zero orbital 
diamagnetism. In this connection it will be apt to clarify 
the case of a circular ring which is also an unbounded sys- 
tem, and yet the associated orbital diamgnetism is readily 
seen to vanish 12 . Is this a counter-example? The answer 
is a definite no. The vanishing of the orbital diamagnetic 
moment in this case is due simply to the fact that here the 
only allowed motion is necessarily along the tangent at any 
point on the ring. The Lorentz force, therefore, must act 
in a plane normal to this tangent. The latter can have no 
effect whatsoever on the tangential velocity because of the 
holonomic radial (rigid) constraint that balances all radial 
forces. Thus, the magnetic field generates no orbital re- 
sponse for this constrained system. This point has also 
been noted by [?]. Clearly, for any non-zero effect of the 
magnetic field on the orbital motion at all, the tangent 
space must be at least two dimensional, as for the case of 
a sphere for example. Thus, for the case of a ring (or for 
that matter, for any non-planar loop in general) the mag- 
netic field is irrelevant, and these cases do not constitute 
valid counter examples. 

The surface of a sphere is, of course, rather special in 
that it has a constant curvature, and the avoided can- 
cellation may well be a consequence of this, it may be 
argued. This has further motivated us to consider in the 



present work the more general case of a closed surface of 
non-constant curvature, namely, an ellipsoid of revolution. 
We consider a charged particle of charge — |e|, mass m 
moving on the surface of an ellipsoid of revolution around 
the z— axis with semi-major axis a and semi-minor axis 
h = at, where e is the ellipticity. The position of the 
particle can be described by the polar angle 9 and the az- 
imuthal angle (/). The Langevin equations of the particle 
in the tangential (0— like) and the azimuthal direction are 
given by 

(e^ sin^ 9 + cos2 9)9 + sin 9 cos 9 [(e^ -1)9'^- 02] 

V sin^ 9 + cos^ 9 
ujr sin 9 cos 9 ■ 



7 sin^ + cos^ 9 



Ve2 sin^ 9 + cos^ 99 + y^] fe{t), 



(17) 
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2 cos 9 9 6 = 



= — cos W f — sm 
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where 77 — 2kBT / (ma'^j'^). The notations above are again 
the same as in our earlier paper [5] . The area normalized 
magnetic moment 



(sin^ 9 ( 
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(19) 



As is clear from Fig. 1, it gave, to our pleasant surprise, 
a diamagnetic moment that now depended systematically 
on the ellipticity parameter e. This case too confirms the 
avoided cancellation responsible for our orbital diamag- 
netism. 

In conclusion, the classical Langevin dynamics of a 
charged particle moving in a magnetic field on an un- 
bounded surface is found to give non-zero orbital diamag- 
netism as found in our earlier paper [5]. We have suc- 
cessfully provided a rebuttal to the objection raised in [t] 
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Further, we have generalized our earlier results 
to surfaces of nonconstant curvature, namely an ellipsoid 
of revolution. The results clearly suppport the idea of 
avoided cancellation. The diamagnetic moment is found 
to depend systematically on the ellipticity. 
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Fig. 1: (Color online) Plot of the area normalized magnetic 
moment, /le eqn. ( |19[ ), vs. the ellipticity e for an ellipsoid of 
revolution. Here uJc/'y = ±5.0 and r] = 1.0. 
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